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We study finite temperature phase transition of neutral scalar field on a fuzzy sphere using Monte 
Carlo simulations. We work with the zero mode in the temporal directions, while the effects of the 
higher modes are taken care by the temperature dependence of r. In the numerical calculations 
we use "pseudo-heatbath" method which reduces the auto-correlation considerably. Our results 
^s^j , agree with the conventional calculations. We report some new results which show the presence of 

meta-stable states and also suggest that for suitable choice of parameters the symmetry breaking 
transition is of first order. 
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I. INTRODUCTION 



O-T n n 

QFT's on non-commutative spaces have been studied from various perspectives recently [1|, 2, 

S 0, IE B S II] • Most frequently studied NC space is the well known Groenwald Moyal space j| 
R 2 e d and various issues like, renormalisation, causality, solitons, statistics have been analysed in the 
literature [l(| EH E EU • The conventional quantisation of fields on these spaces have led to an 
interesting behavior known as IR/UV mixing. The phase structure of fields on such a space reveals 
I/"") ' a new phase known as strip phase [l4j . Alternative quantisation which preserves a twisted Poincare 

0^ ■ symmetry in these theories avoids such a difficulty [IH, [l6| . 

On the other hand the fields on fuzzy spaces like fuzzy spheres, fuzzy CP n etc are explicitly finite 
and do not have the fR/UV mixing U3|. But there is an anomaly in the finite case which reveals 
\£) ■ itself as generating the IR/UV mixing. There is lot of confusion about taking the limit of continuum 

in these models and it has been pointed out various possibilities do exist [13, [H, [H, H3| ■ 

The QFT on fuzzy sphere is a matrix model and easily amenable to simulations and numerical 
studies [2ll I22I [23(1 . We study a real scalar field on the fuzzy sphere using Monte Carlo simulations. 
The earlier studies involved metropolis algorithm to ensure the randomness of fluctuations but 
the autocorrelations are reduced using over-relaxations [23[. But we will use another technique 
extensively used in the study of Higgs model - known as pseudo-heatbath [24]] . Using this algorithm 
we have been able to reproduce previous results from different studies. Apart from this we are 
able to characterise the order of the transitions between the different phases. In particular we find 
the transition between order <-> non-uniform transition is of first order which is mostly due to the 
presence of many meta-stable states at low temperatures in the model. Also we find new results for 
the structure of the phase diagram as well as for scaling of the location of the triple points in the 
continuum limit. 

The paper is organised as follows: the sec(2) introduces the QFT on fuzzy spheres; sec(3) discusses 
the pseudo-heatbath technique and brings out its salient features. Sec(4) reproduces the known 
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results and discusses our results on the nature of phases and the transitions. In sec(5) we present 
our conclusions. 



II. THE MODEL 

We use the following action for the massive, neutral, scalar field on the fuzzy sphere of radius R 

PIED, 11, 

Air 

S{$) = — Tr [<& [Lt, [Li, *]] + R 2 (r<f> 2 + A$ 4 )] . (1) 

Here $ £ Mat^v is a N x N hermitian matrix. The first term in the action is (kinetic) coming from 
the variation of the <E> field on the fuzzy sphere. The quartic term represents the self interaction of 
the $ field. For the thermal behavior of the $ field one needs to study the system in 2 + Id. However 
one can consider the above action as the dimensionally reduced version of a 2 + 1 dimensional action 
with the effects of temperature going into the temperature dependence of r. Finite temperature 
behavior of the <& field with fluctuations included is then studied at different r. In the mean field 
approximation the expectation value of $ field which minimises the action is given by, 



(*) = ±01, = 




<p > for negative values of r and zero for r > 0. The system is in ordered phase for r < and in 
the disorder phase for r > 0. So at the mean-field level there are only two phases. ($) decreases 
continuously to zero in the limit r — > 0. At r = the system undergoes a second order phase 
transition with mean field critical exponents, (3 — 1/2, a = etc.. In the disorder phase (<£>) = so 
the Z2 symmetry of the model is restored. Even though the above form of ($} minimises the action 
there are additional local minima or meta-stable states. For these states, form of ($) is non-identity 
in general. Some of these states will have same potential as the ground state. In the mean-field 
approach these states do not play any role in the phase transition. However they become important 
when fluctuations beyond mean-field are considered. 

Hence the next step in the calculations is to consider effect of thermal fluctuations beyond mean 
field. It is important then to ask if the results of the mean field analysis survive. One expects that 
the fluctuations will destroy a non-zero ($) even for r less than the mean-field critical value which is 
zero. Further more the fluctuations may lead to new phases and different types of phase transitions. 
These are some of the issues of intense numerical investigations recently [19l |2~lI . [23j . So far the 
results show the appearance of a new phase called non-uniform ordered phase. These studies are 
mostly done using a standard Monte Carlo simulations with metropolis algorithm. 

In the present work we study the fluctuations in the above model using a different numerical 
technique known as "pseudo-heatbath" method. Like the previous studies we also observe the non- 
uniform phase. However our results seem to indicate that there arc phase structure within the 
non-uniform phase. These phases can be probed using different operators/order parameters. As a 
consequence there will be multiple triple points in the A — r plane. In the following we describe the 
"pseudo-heatbath" method. Subsequently we present and discuss our results next section. 

III. NUMERICAL TECHNIQUE 



Effects of the fluctuations beyond mean field are computed from the partition function, which in 
the path integral approach is given by, 
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Z oc / D<5>e~ s ^\ (3) 



The standard numerical methods adopted for this integration are Monte Carlo simulations. In the 
Monte Carlo algorithms, one generates an "almost" random sequence of $ matrices by successively 
updating elements of <!> taking into account the measure and the exponential in the integral above. 
This sequence of $ is then used as an ensemble for calculating averages of various observables. For 
example, thermal average of $ is given by, 

1 N m 
m— 1 

here, 4> m is the mth element of the ensemble. Usually there are different ways t o g enerate the 
ensemble. Previous studies of this model have considered the metropolis algorithm [2 lL [23| . In the 
metropolis updating usually there is a substantial correlation between <E>'s in the sequence. For a 
good ensemble the auto correlation between the configurations in the sequence must be really small. 
Though this auto correlation can be reduced by using some over relaxation programme [23]. The 
auto-correlation is greatly reduced, however, when "heatbath/pseudo-heatbath" type of algorithms 
are used [24]]. This method is very much common in the non-perturbative study of $ 4 theories in 
conventional lattice simulations. It gives better sampling and is efficient at least for smaller A values. 
This is why we make use of "pseudo-heatbath" technique. In the following we explain the algorithm 
in greater detail. 

In the "pseudo-heatbath" algorithm, given a $ we update the elements of this matrix one at a 
time. Advantages of updating matrix elements were demonstrated in Ref. [25]. Keeping in mind 
that $ is hermitian we update 4>jj and <F,i simultaneously. We update $y using the probability 
distribution, 

where = a - A) 2 + XB ($y — C) 4 , (5) 

A, B, C may depend on the elements of $ (other than a is a parameter chosen so as to 

maximise the efficiency of updating. In the first step a random number is generated using the 
distribution, 

In the second step the newly generated random number is accepted or rejected using the second 
term of S In our calculations we get for some choice of the parameters, in particular small A 

acceptance rate up to 95%. Over relaxation can also be easily incorporated into this algorithm. In 
the over relaxation process we flip the element <f> — $y in the following way, 

cj)' = A - 20 (7) 

then accept it with the probability exp(— SS). SS is the change in action due to flipping. For small 
A this amounts to changing $ by large amount with only a small change in the total action. Even 
without using the over relaxation method we get very small auto correlation. In the following we 
present and discuss our numerical results. 
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IV. NUMERICAL RESULTS AND DISCUSSION 



To study the phase diagram and transitions we make measurements of various observables such 



as 



Tr($), Tr($ 2 ), Tr(S) 



(8) 



at various values rR 2 for different choices of (N, XR 2 ). In order to check our algorithm we considered 
some of the parameters used in previous calculations pjj [2l|, [H| , and found that our results match 
reasonably well with previous results. The results also agreed with mean-field away from the tran- 
sition point. In Fig. 1(a) we show Monte Carlo history of Tr($) for N = 2, R = 1.0, A = 0.63662 
and r = —1.530502 as used in (23|. Without using the over relaxation we get the quality of data 
similar to that of [23| shown in Fig. 1(b). The average values agree but we observe larger fluctuations 
(Fig. la). 
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(b) Monte Carlo history of Tr($) from 
FIG. 1: Comparison between pseudo-heatbath method and Metropolis method 



Order <-» non-uniform transition 

Having reproduced some of the previous results we considered various values of the parameters 
(N 7 XR 2 ) to study the phase diagram. For large values of XR 2 we found multiple transitions [T^ |. 
For low temperatures, or r -C 0, the average value of <I> is essentially identity matrix indicating 
the ordered phase. For larger values of r the average of Tr($) vanishes while the average of some 
elements <t>y is non-zero. Such a form of <I> average indicates a non-uniform phase which breaks 
spatial rotation spontaneously. Since the average of Tr($) is non-zero in the order phase and zero in 
the non-uniform phase one can use it as an order parameter for the order <-> non-uniform transition. 

As for the order of the non-uniform <-> order transition we find it is strong first order for larger 
XR 2 . This can be seen from the hysteresis effects of Tr(<fr). In Fig. 2(a) we show the hysteresis loop 
of Tr($) for the set of parameters N = 25 and A = 0.8. The value of r corresponding to the middle 
of the hysteresis loop is take to be the critical(transition) value r = r\ for this transition. By doing 
simulations for different values of A we find that the strength of this order <-> non-uniform transition 
varies. For smaller A the transition becomes weaker. For example, for A = 0.4 the hysteresis loop 
was not prominent. 
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(a) Hysteresis loop of Tr($) for TV = 25, R = 10.0 and (b) Monte Carlo history of Tr($) for N = 25, R = 10.0, 
A = 0.8. A = 0.4 and r = -0.75. 



FIG. 2: First-order Phase Transition 



The history of measurement of Tr(<5) and its distribution is shown in Fig. 2b for r close to the 
corresponding critical value n. One clearly sees three degenerate ground states here. Out of these 
two are connected by Zi symmetry and the 3rd has Tr(<f>) peaked around zero and represents a 
non-uniform ordered phase. This distribution indicates that the order and non-uniform phases do 
coexist suggesting first order nature of the transition between these two phases. Even though the 
transition is first order, it is weaker compared to the previous example. 



Non-uniform <-* disorder transition 



In the non-uniform phase Tr($) keeps fluctuating around zero. One can clearly see measured 
values of Tr(<I>) form bands symmetrically situated around zero as shown in figure Fig. 3(a). The 
band structure are not seen in the Monte Carlo history of 7>(<I> 2 ) and S(Q). This band structure of 
Tr($) we saw mostly in the case when there were many meta-stable states before transition in the 
ordered phase, i.e for larger N. In Fig. 3(b) we show the histogram of Tr(<I>) which clearly shows 
a peak close to zero. This implies the state with lowest Tr(<f>) is the ground state of the system 
and other bands are meta-stable states. The meta-stable bands tend towards zero as we increase r 
further. At the same time some bands disappear and/or others merge with the middle one. 

In the basis we choose to work with even though Tr($) fluctuates around zero both $n and $nn 
fluctuate around non-zero Z^ symmetric values. So the symmetry of $ is not restored yet. When 
we increase r the non-zero values around which the first and last diagonal of $ fluctuate approach 
smoothly to zero. Beyond certain value of r = T2 all elements of $ fluctuate around zero restoring 
the Z?, symmetry. In Fig. 4 we show the distribution of these elements of <E> both below and above 
r%. Given this behavior of $11 and $>nn one can consider any of them as the order parameter for 
the non-uniform «-> disorder transition. This also implies that the higher spherical harmonics are 
becoming important for this transition in our basis [19j |. When is approached from below the 
peaks of the distribution of <5>u,$nn smoothly approach zero indicating only one ground state at 
any particular value of r around the transition point. So this transition is a continuous transition. 

Phase diagram and triple points 

When we analysed the data for fluctuations of Tr($ 2 ) we found these peaked at a certain value 
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FIG. 3: Monte Carlo history and histogram for n < r < r<z 
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(a) Distribution of $n for N = 12, R = 10.0 and 
A = 0.04. 



(b) Distribution of ® NN for N = 12, R = 10.0 and 
A = 0.04. 



FIG. 4: Non-uniform «-> disorder transition 



of r within range r% < r < ri (Fig. 5b). This suggests finer structure or phases in the non-uniform 
phase. The peak in Tr($ 2 ) then corresponds to the transition between these phases. 

We did not see any dramatic change in the variables such as Tr(<I>) or the diagonal elements of 
$. The finiteness of the peak implies that the transition is a continuous transition. For larger A the 
transition point was far from both n and Ti. We anticipate that the non-uniform phase has more 
structure than what we see from the behavior of Tr(<f> 2 ). This finer structure could be explored 
by appropriate operators such as multi-trace operators. Note that this fine structure becomes more 
prominent for larger N, which implies that it will survive in the continuum non-commutative limit. 

As r is increased from some large negative value the system explores all these phases for large A. 
For smaller A some of these phases will not appear when r is varied. This leads to presence of triple 
points in the A — r plane. For some small A there is transition directly between order <-> disorder 
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(a) Susceptibility of Tr($) for N = 25, R = 10.0 and 
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(b) Susceptibility of Tr($ 2 ) for N = 25, R = 10.0 and 
A = 0.9. 



FIG. 5: Transition between non-uniform phases 



phases, i.e |ri — r2 1 vanishes. This leads the triple point which has the lowest value of XR 2 . In Fig. 6a 
we show the phase diagram for N = 25 in the XR 2 vs rR 2 plane. This is the triple point studied 
in previous works. In these studies the triple point was obtained by using numerical results for the 
order-non-uniform transition and the analytic results which takes into account only the potential 
term [l9j . In our case both transitions lines are from our simulations. 

Conventional lattice regularization of the model does not show any evidence of non-uniform or- 
dered phase. So it is imperative to study what happens to the non-uniform ordered phase in the 
continuum limit. If the non-uniform phase survives this limit then only it can be physically relevant. 
We have studied the N dependence of the triple point. Since it's not practical to do simulations for 
very large N, one must study scaling to find out the limiting position of the triple point for larger 
N values. In the Fig. 6b we show value of Y tri = XR 2 , X tri = rR 2 corresponding to the triple point 
for different N. The values of N considered for our simulations are N = 4, 8, 12, 16, 25. When 
N was increased the triple point moved away from the origin. Our results suggest that the triple 
point scales with (N 11 , N v ) with /i ~ 1.0 ~ v. We did not observe any universal scaling of the phase 
boundaries in the phase diagram for different N. 

For smaller XR 2 there is only one transition, the order <-> disorder transition. For larger N 
the distribution of the observables such as Tr($), $n, &nn, close to the critical point, show a 
plateau around zero with highly non-gaussian features. This can be seen in Fig. 7 where have 
shown the histogram of TV($),$n around the critical temperature. The parameters considered 
here are A = 0.4, R 2 — 10 2 , N = 12. We take the plateau structure around zero as an indication 
of a transition which is stronger than second order transition. However it does not rule out the 
possibility of second order phase transition for smaller values of A and N. 

Meta-stable states 

We also studied cases with very large values of XR 2 . When r is large negative we find different 
average values for Tr(<I>) for different initial choices of <&. These different values correspond to local 
and global minima of the effective action. The barrier between these states inhibits the transition 
amongst them. The number of these states which we observed grow with N. This can be more 
or less seen from the analysis of the action itself as fluctuations are not much important for small 
r. The state with highest Tr(<E>) found to satisfy (<&) oc 11 and has the lowest value for the action, 
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FIG. 6: Phase diagram and scaling of the triple point 
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hence is the ground state. So we conclude that the state with largest Tr(<5) is the global minimum 
of the system. Other states, which are basically the non-uniform phases, are meta-stable. We think 
that the meta-stability increases with decrease in the average of Tr($). In Fig. 8(a) we show a brief 
Monte Carlo history of after thermalization, for different initial <i>'s. For higher r the bands persist 
but move slowly towards zero. After certain value of critical r most of these states are observed only 
for sometime in the Monte Carlo history and then the values of Tr(<i>) jumped to zero as seen in 
Fig. 8(b). 



V. CONCLUSIONS 



We have developed a "pseudo-heatbath" algorithm to study the finite temperature phase tran- 
sitions of <i> 4 theory on a fuzzy sphere. The results from Monte Carlo simulations clearly show 
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(b) Monte Carlo history of Tr(0) at iV 
A = 0.7 and r = -0.885. 
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FIG. 8: Monte Carlo history for r < n and for r ~ ri 



finite temperature transitions. For some range of Ai? 2 , in particular, for large values one clearly 
sees stable non-uniform phases for some intermediate temperature, intermediate values of r. The 
various phases are characterised by different properties of $. In the ordered phase this behaves like 
a identity matrix. All non- uniform phases have zero TV($). Their existence is confirmed by the 
peak in the fluctuation of Tr($ 2 ). Tr($) serves as an order parameter for the order <-> non- uniform 
transition while <j>n and &nn describe the non-uniform «-> disorder transition. 

The order-non-uniform transition is found to be first order. This transition was found to be strong 
first order for larger values of XR 2 . We conjecture that the first order nature of the transition has 
to do with the presence of meta-stable states discussed above. In fact the state with 7>(<I>) = is 
meta-stable for small temperatures when $ oc 1 is the absolute ground state. Fluctuations can only 
stabilise if it is a stable configuration at higher temperature, so there is always a barrier with the 
ordered phase, leading to first order transition. For smaller values of XR 2 the barrier between the 
stable and meta-stable phases is not high so thermal fluctuations make $ hop between the different 
states. In this case we rather study the distribution of Tr($) to infer the transition value of r for the 
order <-> non-uniform transition. From the distribution we find that ground state is discontinuously 
changing. Moreover the distribution of Tr(<I>) is very non-gaussian, rather flat near zero, suggesting 
that there are degenerate states characterised by zero and non-zero values Tr($). So from our 
results the transition of ordered phase to other phases is always first order for the parameter space 
we have explored. All other transition in model appeared to be continuous transitions. 

The results from previous studies have shown that by doing simple scaling phase boundaries for 
different N coincide [l9| . This expected scaling with N does not occur up to the value of N we have 
studied, though it is the largest so far. We studied the behavior of the triple point for large N and it 
scales approximately linearly with N. However our results seem to agree with the previous studies 
in that the non-uniform phase survives the continuum limit. In our analysis we considered primarily 
Tr($), Tr($ 2 ), $n etc.. However analysis of the full matrix may result in better understanding of 
the phase structure, such as variants of the non-uniform ordered phases. 
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